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ABSgRAOQ} 

A picture can be derived from its "Higihs” 
and "Lows”. "Lows" need very little bandwidth for 
transmission. Hence an efficient coding of "Highs" 
will ensure a high redundancy reduction in the 
entire signal domain, This can be achieved by 
synthesising the "Highs" from the contour infor- 
mation. 

A novel method of coding has been suggested 
for such contour information using Spline Approxi- 
mation . Application of this' encoding procedure 
to "Synthetic Highs System" results, in a large 
compression in bandwidth required for image trans- 
missicn. This method of encoding the contour may 
be directly extended to contour maps, smooth line 
drawings and similar figures. 



Chapter I 
IIWRODUCIIOT 

1 .1 MOO:i7A!DI01T FOR PICTURE RATA COMPEBSSIOR 

There has been a ever-growing demand for trans- 
mission of pictorial information in recent years, 
particularly after the advent of satellites* Introduction 
of satellites in the scientific world has opened up 
many new avenues for exploration and revolutionised many 
areas of application. Consequently, there has been an 
explosion in the volume of data flow all over the world. 
Even such diverse fields as mass media communications, 
space exploration, these days, require enormous amount 
of information, mostly in the pictorial form, to be 
transmitted over large distances* Hence, much attention 
has been given, in recent years to the field of data 
compression, in particular pictorial data compression. 

The reasons for the above is two-fold. First and 
foremost is the heavy flow of information, to be handled. 
Besides the problems -due to voluminous amount of data, 
the large distance to be covered imposes stringent 
restrictions on the quality of transmission by conven- ; , 
tional systems . Digital systems like PGM do provide a , 
high quality tranSmisSipn, 'unaffected by distahce. But , 
the bandwidth re quiremen ts of such a sys tern is unusually 



high. Hence, considerable interest has been shown in the 
last few years towards developing a highly efficient method 
of coding the image information. 

1.2 HACUEE OP PEBVIOUS STUDIES 

The investigations in the field of bandwidth 
compression could be classified into two categories: 

(i) Investigation of highly complicated and 
sophisticated systems with a view to obtain a very large 
reduction in transmission rate. 

(ii) Study of simple practical techniques giving 
moderate bit reduction. 

The reason for studyii3g highly complicated systems 
is to obtain a sort of an upper bound on bandwidth reduc- ; 
t ion ratio that can be achieved without appreciable 
degradation in subjective picture quality. Sven though 
some of these may not be implemented easily at the present 
level of 'art*, such limitations were not taken into 
account in view of the objectives specified. Because of 
the complexity of such systems, these have been studied by 
means of computer simulation. Purther, a computer Simula-; 
tion has the great advantage of flexibility and consequently 
permits an extensive analysis of the entire System. 
of the studies done so far has been limited to monochrome 
still pic tuts s since simulation of motion picture requires 
a large storage capacity. The yarious extensions of 



Synthetic Highs to two dimensions fall under this categoiy. 

The Delta modulation scheme and its modified versions 
like st modulation constitute the class of simple, practical 
systems. Although these have been highly successful for 
speech transmission, these have not produced satisfabtoxy 
results in transmitting picture signals. 

1 .3 BRIEP EE7IEW 

The thesis under discussion falls under the 
category of complicated systems mentioned above. Most 
of the work on this subject, has been done by the 

7 15 

Cognitive Information Processing Group at M.I.T. ’ . 

Since the system under study is closely related to 

1 1 7 

Schreiber's Synthetic Highs System ’ the following 
discussion will be confined to papers related to above 
system only, 

1*3»1 Bandwidth Compression of Picture Signals: 

•Basically, the reduction in bandwidth is possible 
due to the following factors: 

(i) The high bandwidth allowed for conventional 
TY pictures is mainly needed to reproduce the occasional 
shayp changes in intensity, occurring near the contours 
Since the , contours occupy relatively small portion of thS ;• 
picture, there is no need for providing such large band-, 
width for transmitting information whose variation is ’ 



slow except at relatively small number of points^ i'*e., 
statistical properties of the pictiire should be exploited. 

(ii) Human eye does tend to distort the picture 
in its own way. Hence any attempt to distort the picture, 
without lowering the subjective quality of the picture 
is welcome , i.e., the visual properties of the eye are 
also to be taken into account. 

Schreiber’ s philosophy is to treat the picture as 
a dfunction of several dimensions, .spatial as well as 
temporal. It has been found that such an analysis can 
result in a much higher reduction than the ons obtained 
by processing the video signal output of the scanner, 

This is due to the fact that multi-dimensional treatmeiit 
of the picture is capable of exploiting the interframe 
and interline correlations, 

1 *3 *2 Synthetic High Systems 

Experimehts indicate that the human eye tends to 
emphasise the edges in a picture but is relatively in- 
Sehsitlve to amount of changes in brightness over edges. 

On the other hand, in areas where brightness changes ' 
slowly, qu^tization noise is easily discernible . inhere- 
fore, the edges and slowly varying part of the picture 
are to be treated differently , Figure 1 .1 shows a block 
diagram of the system, Ihe video signal, derived from 
the picture by scanning, is pas sod through a low pass 
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filter with frequency response L(^'w) . If the bandwidth 
of the low pass filter is 1/1 0th of total bandwidth, 
the sampling rate has also to be reduced by a factor of 
TO. The signal S(x) is passed through a differentiator. 
The output ^ contains mainly edge information. If the 
low pass signal a(x) and the edge ^ are transmitted 
exactly, the original signal S(x) can be synthesised by 
passing through a filter of response 

H(3w) = (1 - I(3w))/3W (1.1) 

and, adding to it, the low pass signal a(x). This system 

Q 

has been extended to tv/o dimensions by Pan and later on 

6 

improved further by Graham . 

In Pan’s extension of Schreiber’s system, the 
edges are determined by a local threshold function and 
information about the position of edge points v/ero 
transmitted after approximating these by straight lines. 
No local operator similar to the derivative in single 
dimension was used. A reduction of 10:1 to 20:1 Viras 
obtained. Pictures obtained, however, were of rather 
poor quality, although still looked reasonable, 

'in 

Schreiber suggested a mathematical extension 
of his scheme to two dimensions. The b^^^ 
such a scheme is given ; in .Pigute 1 .2 . In this system > 
two different local operators were suggested for edge 
detection. They are (i) Gradient and (li) iapiacian. 











This systom was later simulated by Graham.- It was soon 
apparent that gr-adient operator is better suited for edge 
detection, since it is less susceptible to errors due to 
thresholding than Lriplacian, The pictures thus obtained 
were of extremely good quality with reduction ratios 
varying from 4:1to23!l. 

The picture ern be obtained directly from the 
gradient alone. However, in such a system even a small 
error in re construction will generate highly objection- 
able noise patterns, sometimes masking the picture ■ 
itself. The presence of IPF reduces such effects (see 
Appendix I) due to clnmnel noise and other errors,.’. 

1.3.3 The Proposed System! 

The system under study combines the approximation 
techniques of Pan and the gradient reconstruction proce- 
dure of Graham. However, appreximation is done using 
cubic splines instead of stroi^t lines. This not only 
gives a substantial reduction in bandwidth but, in 
addition, Smoothens the contours as v/ell. The system has 
been simulated on IBM 7044 and reduction ratios ranging 
from 6:1 to 40:1 can be obtained. 

1*4 HJBSIS OUTIIIE 

The Second Chapter cbmprises of an introductioh 
to spline functions aloiog; with a discussion of various 
algorithms of spline approximation, , In the third Ohapfer-, 



the salient features of the proposed systems have been 
discussed in a detailed naoner. She fourth Chapter deals 
T/ith various aspects of simulation. In the final chapter 
a critical evaluation of the system performance is 
presented. 



Chapter II 
SPLIHE APPROZHailON 

2.1 WHAT IS A SPLIIIE ? 

Ever since the days of line drawings, nan had to 
tackle the problem of fairing a smooth curve through a 
sot of points. Often, this v/as accomplished through the 
use of long, thin and elastic strips of wood, anchored 
in place by attaching honyy ’’ducks" at points along the 
curve. Smooth lines, so obtained, are the forerunners 
of what are now widely known as "Spl-ine Punctions" , Thus, 
a Spline is just a smooth curve or rather an assembly of 
pieces so joined as to present the smoothest possible 
appearance* It is this smooth disjointed nature, more thm 
anything else, which makes it an automatic choice for 
approximating contours. 

2 .2 SPIiIKBS AS APPEOXIllATIlTG PURCTIOIS 

It is often acknowledged that pol 3 mcffiials and 
similar ma themat ic al f uncti ons are inade quat e to approxi- 
mate functions from the physical world* Moot of the time , 
the physical relationships could be more effebtively v 
represented by ftmet ion of a disjointed or disassociated - 
nature than by the more familiar functions of mathematical 
world. One of the chief advantages of a Spline is that, its 
behaviour in one .section can not determine its subsequent 



course in another section. Besides , the "Enots" 
enjoining the various pieces, can be varied in position, 
thus providing a wide variety of smooth curves . In fact , . 
much of the interest in splines arises from our ability • 
to choose or vary its knots. Above all, Spline has an 
inherent smoothness associated with it as a natural 
consequence of its minimum curvature property-. 

2,3 kAIHEMAIICAL BSPElSBlffiATIOlTS OP SPIIHBS 

A precise mathematical description of spline is to 
portray it as a piecewise polynomials of degree n that aie 
so connected as to have (n-1 ) continuous derivatives at the 
joints, i.e., it belongs to class of fimctions denoted by 

Spline functions aro usually defined in terms of 
their knot set given by 

0»kQ<k^<. .. k^ < k^_^^ *1 \ 

The representation of spline function is given as follo7/s: 

S(X,E,z) = a.(x-k*)J.+ P(x) (2.1) 

where P(x) is a polynomial of degree n v/ith coefficients 
a^, i != n+1 , m+n+1 and X, K, . (z-kjj_)^ are given by 

55 '■■'1 HI' 

= .'O-'. .otherwise^' 



A spline, v/itli m interior knots, has m+l polynomials 
chained to form a single entity. Since a ploynomial requires 
only (n+1 ) parameters for unique definition, it is possible 
to represent the spline uniquely by (n+1)(m+1) parameters. 
However, it seems from equn*(2,1) that a spline can be 
represented using only (n+2m+1 ) parameters. Ihe additional 
m(n-1 ) parameters normally needed to specify (m+1) poly- 
nomials sections, have been rendered superfluous by -the 
m(n-1) continuity constraints imposed at the joints to 
ensure the smoothness of the spline, Thus, the continuity 
constraints not only introduce smoothness in the curve but 
reduce the number of defining parameters too. 

Ihere are many other fforms of representation apart 
from equn,(2,1). But, the most useful one from the point 
of vieyy of approximation is given; belov/. In this represen- 
tation the position of the knots, value of the function at 
these knots plus the. (n-1 )/2 end conditions at each of the 
ends are used to specify the spline of, an odd degree, That 
is, S(A,K,x) is represented by 


v^^o for i = 0,1 , ...,m+1 
^ 0,3 ^'^ ’■^m+.1,j 3 = 1,2,..,,(n-1)/2 

where 



^ s(r,i,x) 

dx^ 



( 2 . 2 ) 


The above representatioh, by providing direct 


information about loiots greatly simplifies Goding proeedurc. 



Also, it requires no more than the minimum possible 
(n+2m+1 ) parameters. All these factors render this 
particular representation .highly suitable for coding 
purposes, Hov/ever, it has one drawback, viz, the 
reconstruction of the spline becomes laborious. 

To simplify the reconstruction, an intermediate 
representation which can be obtained from equn,(2.2) is 
made use of. In this representation, S(A,K,x) is 
represented 'bj the knot set K and the matrix 

{ Vi^ , i = 0,1 , . . .,m+1 

2 = 1,2,...,(h--1)/2 

where 

^ S(X,K,x)f (2.3) 

dx** l^=^i 

The representation given in, equn,(2,3) gives all 
the necessary information about each segment of the 
spline. This permits each segment of the spline to be 
treated indepenaent of others, thus lessening the 
computational effort involved in the reconstruction of 
the spline , It is easy to see that any such computation 
without arriving at the parameters specif ied in equn.(2 ,3) 
directly from equn,(2,2) is difficult and laborious; Th^ 
representation given in equn. (2 .3) .can be easily obtained 
from equn. (2 ,2 ) by solving the vappropriate^ m(n-i) conti- 
nuity equations. In the case of cubic splines, this 




reduces to the problem of solving a system of simultaneous 
equations, represented by a tri-diagonal matrix. 

2 ,4 SPIin-IE APPROXI]iIAO?I01I ilESEHODS 

There are many algorithms of spline approximation 

5 12 

available in the literature ’ . These can be broadly 

classified into two classes. They are; 
i) Approximation with fixed knots, 
ii) Approximation with free or variable knots. 

Although the second one seems to be of relevance 
to ttie present problem, the first one will also be consi- 
dered in detail, since all methods concerning the second 
case need a fixed knot approximation technique as a 
subroutine , • 

2.4.1 Pixed Knot Approximation Problem: 

This class of approximation problem may be 
subdivided into two subclasses based on the norms used 
for approximation errors* Two of the most popular ones 
used are the norm or the minimax error criterion and 

the L2 least squares error criterion, 

■ ■ 1 g 

■ Schumaker has given a computational method for 
approximation with norm based on exchange process of ; 
Eemez y^iety. Such an approach which assumes an error 
curve qf ” normal" fom with h+m+1 alternati]^ extrema 
suffers a disadvantagje because of the non-Chebyshevian : 
nature of the approximating function. Por some functions, 



the "best approximations have fewer than n+m+1 alternative 
extrema* Perhaps a greater drawback is observed in the 
form of additional error oscillations when joints are near 
their optimal locations. 

2 

Barrodale and Yoiing have used a linear programming 
approach, which thou^ closely related to exchange processes 
of Heme z variety, is not based explicitly on any characteri*" 
sation process and avoids these difficulties* 

J .R. Rice and Oarl Be Boor^ have evolved a proce^- 
dure for approximation which is much more elegant than 
ai3y of the methods described above. He has treated the 
problem as a linear approximtion problem and solved it 
using an. orthogonalization procedure. 

11 8 

Powell and, lawson have also dealt with the above 
problem and a concise discussion of various methods has 
been presented by Esch and Eastman in their paper^. 

Out of the methods described, Rice’s procedure is 
extremely economical in time * Because of the orthogona- 
lizatioh procedure it can be so manipulated that variation 
of a particular knot will affect only one orthogonal 
function and hence , the new error , c ould be obtained from' 
the error at the previous stage of approximation* In all 
the other met hods » the: whole procedure has to be: c ompletsly 
repeated even for a small variation of a particular knot. 



Tills in-built facility is particularly useful in saving 
computer time v/hen used as a basis for free knot approxi- 
mation. Because of this advantage, the algorithm suggested 
by Rice has been preferred to others for use in the proposed 
system. 

2 .4 .2 Fixed Knot Approximation by Orthogonalization: 

The problem is to determine such that 

/ (f{x) - S(A*, K, x))'^ dx (2.4) 

Q , . , 

is minimised over the space of vectors A. 

This can be treated as a linear approximation 

problem and there exists a solution for such problems of 

approximation using or-chogonalization procedure. In this 

method, the problem is reduced to one of finding a set of 

orthonormal functions which will form a basis for the gi-ven 

class of approximating functions. It has been shown by 
12 ■ 

Rice that such a basis exists for spline functions. If 
we denote the set of orthonormal splines by ' 

j = 0,1,...,n+m 

For i-g: n 0^ (l!,x ) « (x) 

where is the; legendre polynomial of degree J 

in (0,^-'1;),.- 

For j > n (K,x) , 3 = h+T , . . . n+m 

are'- given- by; the/'-eqhatiohs'-- 



( 2 . 5 ) 


1 



Thus- starting from Legendre polynomials the 'vtiole 
set of orthonormal splines can he determined, using Gram- 
Schmidt procedure . The optimal approximating spline is 
given hy 

_ _ n+m _ 

S(A*, K, x) = x: a^ 0.(K, x) (2.6) 

3=0 3 D 

1 

v/here a* = f f(x) 0^(K, x) dx 
3 0 3 

2.;4-.3 Variable Knot Approximation Problem: 

This is basically a nonlinear approximation 
problem* Here the problem is to determine the position 
of the knots so as to minimise the error criterion. Once 
a fixed knot routine, giving the requisite approximation 
along' With the erron for a given set of knots., is available, 
the problem can be reduced to mini m ization of a nonlinear 
function which maps any knot set to the corresponding 
approximation error, 

*1 Q 5 /I 

Schumaker , Esch^ and Eice^ have all tackled this 
problem in their own ways. But the discussion: here is : / 
restricted to Elce's algorithm due to reasons 'already ' 
mentioned In 

The problem : is to determine:, an optima 1 , X* such , 
that for any given f,(?c) on (a, b) and4>0, 



^ _ _ p 

J (f(x) - S(A*, K*, 3:)r dx (2.7) 

a 

is minimized with minimum possible number of knots. In 
other words, ^ as wc:!! as A* should be determined 
completely, their dimensions as well as their values ^ 

It has been found theoretically that a solution 
to the above problem is feasible. In support of this, 
wc have the following theorem (Ref .12, pp.145)i- 


Set j^(f) = Pf(x)-S(A-^,r*,x) 


( 2 . 8 ).. 


where the spline S(A'^' x) is a best I2 approximation 
to f(x) with degree n and m variable knots, then either 


<2.9) 

Thus, in theory at least, it is possible to 
achieve the approximation to ary desired degree of . 
accuracy with sufficient numbers of knots. However, such 
an approximation with minimum number of knots is of parti- 
cular r el evahee here. This is what Rice and De Boor have 
attempted to achieve bj their eomputational algorithm. . 

2.4;.4 ..Ricel a, ■Method.. . 

J.E.Rice and Carl Be Borr^ have suggested a procedure ' 
for cubic spline approximation with minimum number of free 
knots. This procedure requires as a basis a set of 



orthonoimal functions to be determined* I'irst-n+l 
(here, 4) noimialized Legendre polynomials are obtained 
and f(x) is replaced by . 

f(x) - X a* 0.(x) (2.10) 

j=o , ^ 

She remaining 0^^^ (K,x), j = 1,,2,...,m are determined 
so that 0^^^(K,x) involves only the y leftmost knots. 

Lhis makes it possible to economize on computations later. 
The error can be determined as 

e(x) =s f(x) “ X a. 0. (K,x) (2.1 1) 

3=0 ^ ^ 

Given an initial set of knots, they arc improved 
cyclically in oitler to minimise the error. She Cycle 
starts with rightmost ioaot and working to the lefty 
each is varied so as to reduce the error as a function 
of this knot. This cjrclic procedure is continued till 
suitable convergence criterion is met. The criterion of 
a particular knot is made using discrete Howton method. 
While varying a single knot k^, a temporary orthogonal 
basis for splines > not involving fc^, is obtained . This 
basis is used along vvith one: rothpgonal basis function 
Oj^(x) involving k^^ Hence, only this last function ; / 
changes during the discrete lewton method and one Can ■ 
replace’' the:-a:(x) '■■by 

■'V'"- 'o ■ -b''',.' y„! : ■■■, "V ■ 

e(x,kj.) » r / f(xX^^^^ dx (2.12) 



However, in tliis procedure care should be taken 
against crowding of knots, aliasing of knots etc. 

Under no.rma:. circumstances if the error criterion 
is not met with the initial number of knots, additional 
knots are introduced, They are included one at a time 

optimized cyclically till the required error criterion 
or some other exit criterion is met. A point is deter- 
mined, at the end of each stage, where error is maximum 
and the new knot is introduced midv/ay between the two knots 
bracketting the point of maximum error. The algorithm may 
have to bo terminated sometimes even before sufficiently 
small error is obtained, o'wing to crowding of knots 
around a particular point. A similar exit is forced when 
an unusually large number of knots are required to approxi- 
mate the given curve . Accuracy of this method, like most 
of the other optimisation methods, is very much sensitive 
to the initial choice of knots. This is the major dis- 
advantage of this method. Ho’wover, this disadvantage is 
common to all the other methods as well. 

Although "tiiis algorithm can be extended to splines 
of higher degree, it is not advantageous to use higher 
degree splines because of the larger number of parameters ; 
required to represent , them. Also , the complexity of the • 
procedure increases enormously even for simple problem with 
fixed knots , Hence ; onlj' cubic splines have been used in the 
.■■system' under 'study,. 



Chapter III. 

PEOPOSED SYSTM 

3.1 lEIHODUCIIOE 

Ihis system combines the salient features of Pan^s 

Q 

and Graham’s schemes. Pan demonstrated that contours can 
be coded efficiently 'hj proper segmentation, folloy/ed by 
straight line approximations of the various segments. How- 
ever, Pan adopted an "ad hoc" procedure for reconstructing 
the picture. Consequently, he succeeded in only producing 
a pictui'G which could not quite meet the requirements, 
dera -.mded of it . 

■ ' 

Graham subsequently developed an elegant but simple 
method of reconstructing the picture from its contour infor- 
mation. This scheme could not only impart the requisite 
quality to the received image but achieve this v/ith ample 
reduction in bandwidth, as well.: Although Graham exploited 
the connectivity of the picture contours, he did not mahe 
use of the natural "smoothness" inherent to the contours of 
the picture. Hpr did he effectively’’ use the correlation 
betv;een the direction of contours and the gradient direction 
in siibse quent coding of the contours , prior to the actual ‘ 
transmission,-"' 

In the proposed systcmi , the Graham’s reconstruction 
procedure has been incorporated with a view to obtain good 
quality pictures. At the same time an effort has been made 



towards exploiting tho "snoothness” of the contours which' 
renders it suitable for coding by approximation. Also, a 
scheme is developed for extracting the part of the gradient 
information from the approximated contours themselves, 

3.2 SAIIBMIE PSAfUEES 

Since the system differs from G-raham's scheme 
mainly in coding of the cox-itour information, only that 
portion will be considered here. 

A simple block diagram (Piguro 3.1 ) illustrates 
the various parts of the system. The first block forms 
the contour tracer which extracts all the gradient 
contours in the picture and sends the position information 
of the contour and the gradient along these curves separately, . 
The position information is then subjected to a pro '-smoothing, 
process which renders it suitable for spline approximation. 

The parameters of approximation are then fed along \iath 
the coarsely quantised gradient information to the encoder 
for proper coding prior to the actual transmission* 

The receiver shown in Pig. 3. 2 reconstructs the 
splines approximating the contours and organizes the 
gradient along these reconstructed curves. The gradient ' 
map of the picture, so obtained, is passed through an 
appropriate high.pass filter to be odded on, at a later 
stage, to the defbeussed image recoived at a much less 
sampling rate, A detailed discussion of each subsystem , 

■ is 'given below.,.. 













3.2.1 Contour Extraction: 

The contours are extracted using the same procedure 
adopted by Graham^. Gradient of the picture signal is 
calculated for every picture clement. A threshold is set 
foi* the starting point of the contour. Once this is 
realized, a search is made among its 8 neighbours 
(ligiu?e 3 . 3 ) using another, threshold. This process is . 
continued till there is no continuity in the contours. 

This way, certain spurious contours due to noise can be 
eliminated by 'leaving out very short ones. However, 

■ "-i 

this is not the only algorithm available, Eosonfeld et.al. 
have suggested a highly effective edge detection technique, 
Eurthcr in the case of digital picturyos the gradient can be 
defined in innumerable ways, Pingle"*^ and Sakai”* ^ have 
used complicated versions of such definitions to detect 
edges in noisy pictures. But the one Graham had used, viz. 
arithmetic difference of intensity values at adjacent points 
is extremely simple a rd easy to implement. Since the study 
has been done on pictures with very low noise, it was not 
necessary to go in for highly complex algorithms. 

This method of contour extraction needs sli^t 
modification before its implementation* Because of the ^ 
digitization of the picture the contours of the picture , ' 
instead of= being perfectily smooth ^ consists of a large ; 
number of SBmtli edges, both vertical and horizontal. 
Therefore j it is logical to represent this contour 





in tcnas of these sr.all edges. However, Graham for the 
purpose of achieving high coding efficiency has chosen the 
picture element as a basic uniu constituting the conto’jrs* 
This method of represent at ion, in some situations, causes 
a smrt of spread in the contours, leading to a mix— up of 
edge points in adjacent contours (Figure 3.4 )• Moreover, 
since the extracted contour has to be approximated, the 
particular representation used, should only provide a thin 
contour instead of high coding efficiency. Hence the 
small edges constituting the contours have been chosen as 
the basic unit. Those edges ar. represented in terns of 
one of their ends as shov/n in Figure 3 . 5 . The particular 
end used in the representation of the contour will here- 
after be referred to as "nodo’’ . 

3.2.2 Approximation of the Contours: 

It may be noted that the contours extracted from 
the picture have to pass through a device which approxi- 
mates the extracted contours by splines tc the required 
degree of accuracy and gives as output, the various 
parameters of apprexinatien. The advantages of using . 
splines lies in its smootlniess and its ability to 
represent the contours in: an ef ficient manner ^ But the , 
spline functions, for that matter any mathematical 
function, can approximate on]^ "single valued functions. : 
On the other hand, closed contours are very likely to 
■ occur quite of ten in , a picture . This necessitates 
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s ome amount of ‘pre-processing of the contours derived from 
the picture. IThe pre-processing should be such as to 
retain the position infernation as well as provide this 
in tho form of single valued functions. This can be 
achieved in one of the following ways? 

(i) Approximation after segmentation! In this 
ipethod, the contour is divided into many segments in such 
a waj’’ that in each segment the x-co-ordinate of the nodes 
constituting the contours are either non-decreasing or 
non-increasing. Then each of these segments is later 
smoothed so that ho two coilsecutivo nodes representing the 
contour have their x-co-ordinates equal. This is achieved 
by. choosing the proper ends v.j and v,^ of horisontao edge 
and the succeeding vertic::,! edge as nodes (2’iguro 3»5) 

or choosing the node lying in the middle of horizontal 
portion to represent the flat portions of the segments 
as shown in Figure 5,6, In case the flat portions occiir 
in either the beginning cr the end of the segment, the 
ap pr opri ate end of su ch f 1 at p ort ions arc cho sen as shown 
in Figure 3,6. 

(ii) Par trie approximation; In tMs method, each 
co-ordinate is expressed as a, function of contour. lengths, 
a monotonically increasing parctiiqter and approximation: 
procedure is applied to the newly obtained functions . This 
method obviously suffers from the disadvantage that each 
curve has to be effectively .approximated tv^ico besides: 



nccdiiTg double the number of parameters for representing 
the approximation. Also, in the case of digital pictures 
vdicrc. contours are represented by vertical and horizontal 
cdgGS, there is no proper definition of length of the 
contour. Perha.ps, this difficulty can be overcome by 
defining the length of rny segment to be proportional to 
the number of snail edges constituting it . But this does 
not solve the- problem either, since it has been observed 
that the reproduction of such approximated contours is not 
quite satisfactory. This is probably due to the digitised 
nature of the me assure chosen for length. Actually in the 
process of the parametric approximation of a digital contours, 
the smoothness of tho contour is not rotained. , thus loosing 
tho very basic property requi'ncd of spline approximated 
functions. Moreover, the spline approximations arc higlily 
effective for contours of medium length. In ordor to 
achieve this efficiency, long contours may have to be any 
v/ay segmented. Thus for very long contours, v/e are not 
ro ally avb idi ng the s e gme nt at ion p roc e dur e even by this 
method. 

In view of the above mo ntionod reasons, the . 
sogmontation approach has boon adopted for tho system ' 

' und or study , 

' 3 . 2 ',. 3 gradient ; Quantizatibnr ■ ■■■' 

The gfadieht information of : the : picture is very 
vital for a satisfactory image reconstruction, fortunately. 



the hunan ejre is not sensitive to the slow intensity 
changes along the contour. Moreover, even at contour 
points, it tends to onphasize only the change in bright- 
ness. In fact, the eyo is relatively insensiti-ve to the 
magnitude of these changes. Consequently one can do 
with a relatively coarse quantization of the gradient 
magnitude. However, while reconstructing the picture* 
its quality may be inpreved, if the gradient is smoothed 
along the contour. For this purpose splines will be 
highly useful, further Ixjcause of the, nature of coding 
of the contours, the gradient dirocticn can be derived 
directly from the crientation of the, corresponding edge 
since the gradient diroction should alw^.-ys be perpendi- 
cular to the edge orientation. 

3 .2 ,4 Receivers 

Ihe receiver part of the system is illustrated 
in Figure 3.2. She functions of its v-:rious subsystems 
have already been briefly discussed in Sec . 3-.2. 

It has been found, based on the results of the 
simulation that the picture reconstructed from coded 
signal is fairly good to look at and has no particularly 
discernible defects. 



OliTDter IV 


SHiULilTIOlT STUDIES 

4.1 SdE ASPECTS CW SniDLATIOU 

The 'entire syston has been studied through 
simulation on IBM 7044. The various aspects relatirg 
to the choice of the different parameters of the system 
will be discussed in this portion of tlis chapter. 

4.1.1 Iiow Pass Eiltcr: 

This forms an integral part of the transmitter. 
This portion has been simulatod using a routine called 
^IiOWPIL*. Given the input imago samples as well as the 
filter response samples, 'IiOWPIIi* merely c?:lculates the 
convolution between the t'wo, i.c., if h(x,y) is the 
filter response and b(x,y) is the .picture signal, then 
the output of VIiOWPIL* is given by 


i(x:»y) = /|b(x-p,y-q) h(p,q)dpdq (4.1) 


In digitized form, this nay be expressed as 


m.“1 m-1 


" S £ ^d“a,k-h ^a»b 


(4.2) 


where m X ia is the size of the filt of response matrix 
and n X hj that of the picture signal ; matrix. 



p 

HowG'ver, one does not need all the (n+n-1) samples'. 
It is enough to trahsnit samples, corresponding to the, 
picture image points , since these n samples are enough to 
I’cproduce the image signal. Other things being equal, the 
choice of the linear filter is mainly dictated by consi- 
derations of computer tine. It is well Imown that it is 
advantageous to use filters v/ith separable kernels, 
particularly, when the conventional convolution methods 
are employed. Further, the above algorithm proves to be 
oven faster than ones using FFI techniques, when the 
filter samples are small in number. With this view in 
mind, a Gaussian function is chosen as the low pass 
filter, i.e. 

l(x,y) = 2^ expC- ) (4.3) 

Although the filter function extends to infinity on 
either side of the origin, the sample values of tliis 
filter fall very rapidly to zero, after first few samples. 
Hence the filter can be truncated after : a finite number cf 
samples without appreciable loss in accuracy. In our case, 
the cut-off point is choson to be 4*^ and the significant 
samples are normalized so that 

, IJlCxiy) dx' dy ■.,.'^:(4'.4)v, 

Where l(x,y) is the low pass filter response . 



(4.5) 


laking the Fourier transform 
l(w^) = exp(-'51_|£- ) 

whore represents the rc.di£il spatial frequency. 

Since more than SOfo of the energy is at Spatial 
frequencies less than where is given by 


W 

rmax 


1.516 


(4.6) 


this itself,: can be chosen as the cut-off frequency. Since 
the sampliiig theorem holds in two dimensions , a sarpling 
r:-tc of 2f will bo sufficient for transmitting 

the "Lows". 


For cr= 4, the sampling rate s= ■ 

That is, .if a 64^ picture is to be transmitted, its 'Lows* 

2 

need only 8 samples. Further, because of the tv7o dimen- 
sional interpolation used, it should be possible to quantize 
it with fewer than 6bits/saEple if the original picture 
is quantized at the same rate. 


The choice of in effect determines the number of 
filter samples as well os the cut-off . f re quonc^^ of the 
f ilt er .. A low cut-off f re quency' would mean a Substantia 1 ■ 
gaih in the number of bits liequircd for transmitting the 
lows , But , this effectively means a higher value of 
and hence a. .larger number of samples are . required to 

represent the filter response, :the3?eby reducing, the 



speed of the routine. This, however, should not be a 
natter of concern, if v/e confine ourselves to sinulation 
of the corresponding analog system. On the other hand, 
if one chooses to study the behaviour of the digital 
version of the proposed system with a view to better 
the existing POH image transmission, this factor has to 
be taken care of. Also, the reduction, so obtained, 
should not be at the cost of any deterioration in picture 
quality, Graham has stated in his paper, that the value 
of O’ = 4 gives satisfactory results. 

4,1*2 Contour Extraction; 

■The algorithm for contour tracing has already been 
explained in the previous chapter. The routine "OOUTOR"; 
achieves the various functions required of the contour 
tracer. One very important parameter to be chosen is 
the threshold sot up for contour selection. Basically, 
the entire scheme is based on the philosophy that the 
high frequency portions of the signeT occurs only at 
limited places of the pictures, namely contours. 

However, due to quantization and other sources of noise, 
a noisy point may be mistaken for a contour point, if the : 
threshold is not sufficiently high. On the other hand, if 
it is made very high, one is likely to miss some details 
altogether. Eurthe^» pait rof the infprraation has already 
been sent by ’lows*. Therefore , the cut-off frequency of 
the filter which determines the information already: 



transmitted should be closely related to the threshold 
for contour selection , which selects the additional 
information necessary for reproduction. A rough analysis 
done in Appendix II supports this view. 

Since the processing has been performed only on 
two level pictures, the choice of the threshold, did not 
pose any problem. It is worthwhile to note that G-raham 
was able- to reconstruct reasonably good quality pictures 
using a contour start threshold of 25 and edge point 
threshold of 6 for 256 level pictures. ' 

4.1.3 Contour Approxii-iation: 

The algoritimi used for the spline approximation 
has been dealt with in detail in the Chapter II. The 
routine ’APPROX' performs this function. However, one 
important aspect of it, which was not covered in 
Chapter II, will bo discussed here. As already explained, 
for a spline reconstruction, ap^t from the information 
about the knots, two more parameters usually in the form 
of end derivatives have to be specified. Since the 
derivatives can assume any value from +cx) to -go , the 
problem quantizing it is not straightforward. One way 
of achieving this is to quantise the angle of inclination, 
which ranges from +90® to "-90® . Hovvevor, that does not 
fully solve the problem,- If a linear quantization is used 
in such a situation, there is a large error in the dorivativc 



when the inclination is close to right angles. This 
v/as observed while studs^'ing the simulated sji^sten on the 
computer and a remedial scheme v/as tried. According to 
this scheme, the range close to 90° is finally quantized 
and the rest of the range is quantized in a coarser 
fashion so that the error is reduced even in extreme 
cases with more or loss the same number of levels as 
before . 

4»1 .4 High Pass Pilter: 

This is used for reconstructing the high 
frequency portion of tlio signal from the gradient map 
of the picture. This forms a part of the receiver. 

Its design is solely determined by the low pass filter 
used at the transmitter (see Appendix III). This part 
of the system is taken care of by the routine ‘HIPIL’ 
which directly computes the requisite response using the 
conventional convolution method. The speed of this 
routine can be improved by adopting the PPT approach, 

4,2 STUDIES OIT COITTGUE APPROXE^ATIOI 

ThO: effectiveness of the approximating rout ine 
has been tested under various conditions , Two different 
sets of data have been used in all these exporimerttSv 
Two weather maps, each consisting of three bontours, 
in the digitised form constitute the: data. The fi^st 
study was done to compare; the results uaing the tw^^ 



approaches of approxinating the contour. The results 
of this study have been given in Tables 1 and 2 and 
figures 4.1 to 4.6. The results clearly indicate that 
the segmented approach is far superior to the parametric 
method both in respect of the number of parameters of 
approximation and the closeness of the approximation. 

A study of quantization of end conditions was 
also done on these tvro naps and the results shown in 
Table 5 and Figures. 4.7 and 4*8 < The results generally 
indicate that a 6 bit quantization gives satisfactory 
results under normal conditions. However, for curves 
requirii:ig steep end conditions, some error is inevitable. 
This tendency 'can be reduced to a great extent by using 
a high resolution data of the sane map. The results 
obtained are listed in Table 3 and shown in pictorial 
form in Figure 4.9. 

After going through all these results, it appears 
that the spline approxiioation is highly effective when 
applied to cuives of length ranging from 10 to 40. For 
shorter curves, usually a proportionately high number 
of knots are required to tie up the various sections. For 
longer ones, : the approximtion procedure not oniy taksB;t 
much time but also ends up with a solution which is 
optimum only in the local sense . The routine takes n 
time ranging from 5 seoa. to 12 secs, for approximating 
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Effect of Quantization of Ena Derivatives 



















a contour, made up of 15 sample points on IBM 7044* “ liiis 
routine is usually terminated when the approsimation error 
is reduced to a value of about 1 .0 in the minimax sense and 
a value of 0,3 in the least square sense, '<7hen the curve is 
unusually kinked, the minimax error is allowed upto the 
value of 2,0 and the least square error upto 0.9. On the 
average, this type of approximation reduces the number of 
points needed to represent the contour by a factor, close 
to 6 j 1 , All figures, concerning the above subject, were 
drawn on the plotter attached to IHd 1800, 

4 .3 EYAIUATICN 01' SYS2SI PBZE'OHJIAIOE 

Preliminary irrvestigations on contour approximation 
favour the incorporation of the spline approximation scheme, 
into the Synthetic Highs System. However, the system per- 
formance Cc?n be truly evaluated only after processing a 
picture signal, • 

4 .3.1 I/O Structure of the Picture Signal; 

ihe conversion' of a picture into tbe digitised 
form accepted by the computer, is often achieved through 
a scanner equipment which converts the picture into the 
desired form and stores it on the tape, I^he- processed 
picture is recorded back oh the tape , ihe eq.uipment is : 
used once again to convert the digital figures into tho 
normal pictorial form * Another f oi-m of computer output 
often accepted when such n facility ia unavailable, is to 
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represent the various intensity levels 1)7 different 
cor.b illations of alphabets. However, such a picture will 
be coarsely quantized. Besides this needs some of the 
lines to be overprinted on ones already printed, Owir^ 
to lack of access to eitlicr of those facilities, only a 
t\7o level picture was processed. Since the main aim is 
to test the suitability of the contour approximation 
tochnique with relevance to picture transmission, it is 
hoped that the results obtained 7/ould not differ signi- 
ficantly from those using a multilevel pictui« . .Even 
though the input signal itself consists of only two levels, 
th;.:; system is assumed to bo capable of transmitting many . 
more levels so that the ” Lows 'V of the picture are received 
'M’ithout significant loss of accuracy. 

4.3.2 Coding the Contour Information: 

Although the spline approximation 
of points to be transmitted, it does so, as the cost of 
removing the neighbourhood relation between adjacent sanplos 
transmitted. 'Whereas in G-raham’s scheme, for every 
additional point to be transmitted only 3 more bits have tO: 
be used, in this scheme a lot more than three is re'quiced.i 
This calls for stringent limits to be imposed on the number 
of bits required to transmit the contours , :This problem 
can be tackled in the foil owing, way: 

In each : contour , only the starting knot is trans- 
mitted with full details, pertaining to its position*. 
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tiiG subsequent knots onlj tiieir relative positions v/ith 
reference to preceding iaiots are transmitted. Moreover, 
since the distance betr/een the adjacent knots will vary 
from contour to contour, a code word indicating the number 
of bits used for specifying these knots could be transmitted. 
This way, some gain can be made on the number of bits 
transmitted. However, this is done at the cost of 
redundancy and therefore sufficient error correcting 
capacity should be imported to the code word indicating 
L'ho munber of bits re presenting the knots. Further, there 
be scric ^short ccntoui’s in the picture, which can not 
be officioxitlj’- approxii-:;.ited . It is advisable to block such 
contours into one group and transmit them separately in the 
normal mminer of contour coding. This will also inipisve 
the coding efficiency of the entire scheme. 

The picture processed by the simulated system is 
shown in Figure 4 .10 and the results of approximation 
tabulated in Table 4 and Figures 4.11 to 4.1 4. The pro- 
cedure for calculating the required nucaber of bits for 
transmitting the given picture is clearly indicated in 
Appendix 17. ' Iho final results indicate that' a . reduction, 
of high frequency informaticn by a factor of about 6 to 40 
is possible with this system, l.e. , it coji provide an 
improvemont by a- factor of' approximately 1 i.7 over Graham’s 
scheme . Sinoe ‘Lows’ nQod\ only a smll fraction of the 
total number of bits used, those figures hold; good even for 
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tho signal in entire frequey domain. Thus the proposed 
syston can provide a reduction in bandwidth by a factor 
var^'ing from 6 to 40 in conporison with 6 bit PGM 
scheme. 



Chapter Y 
COITCLUSIOITS 

This system has ^.efinite inproveinents to its credit 
v/liich can be favourably implemented in Craham’s scheme. C^e 
important gain made in the coding is achieved by generating 
the gradient direction from the contour orientation. This 
sharpens the picture by thinning its contours which have 
already been smoothened by spline approximation, • These 
effects have not been observed on the results obtained, 
probably due to the fact that only low resolution pictures, 
that too quantized to t'wo levels, have been used. 

This system suffers from the drav/back common to 
all Synthetic Highs Systeias. Any picture not only consists 
of contours defined by changes in gray level, but includes 
a few edges characterised by the statistical properties 
of regions lying on either side. Although there are 
algorithms available for extraction of such contours, 
there is no known methods of regenerating signals providing 
the textural properties . Schreiber has indicated as to 
how to attack this problem. The gain in the bandv/idth 
obtained using splines can very v/ell be utilised to' 
accomodate such textural information. . 

Since the systbn laas not yet passed the crucial 
tost of transmitting mult i*-level pictures, any further 
elucidation on the sub yeet is likely to /bs treated as a 
vain boast. As suoli, a fbllow-up using multi-lever 
■ pic tures/ is' 'deemed' ;0.saQhti'al;i,;''-. . 
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Appendix I 

01 PBHB'OMANCB 01 I'Hd RECEITBE TO NOISY INPUT 


The ansaysis is confined to signals of single 
dinciision, since it con be extended to higher dimensions 
o;ithout any difficulty. The receiv..r in Synthetic Highs 
Systen consists of a gradient regenerator followed by a 
!:ii gl'. p a s s f ilt cr . 

let H(v/) be the rosponso of the HU. Then , , 


E(w) = 1— (AIJ) 

where Ij(w) is the f liter rosponoG of the LPl in the trans 
mi t ter . 


let %(w) be the spectral density of noise input 
to the filter and N^-jCw) the corresponding output density * 
Then 


N,,(w) 


i_r ( w ) 


N^(w) 


CAI.2) 


Total noise output power is giiren by 




L- 

2n 


I 


,+cC' 1 fY™\ 2 


1 - Ii ( w ) 


NoCw) dw 


Iiet us assume that 

■■ ■ ;■ ■ -00 ■ W 


'i ,e , is;' white:.' Then 




Nq(w) is unbounded at w = 0 unless 


lin 

w — ^ 0 


1 - I(w’) 
w 


< cr- 


(AI.4) 


The neccss r:;- condition for (AI.4) to held i 


IS 


lin 


w 


. l(w) = 1 


(AI.5) 


If l(t) i.3 the tine x'esp..nse of the LPF 


“hUO ■ 

~ / l(t) exp(- 3 \rb)dt 

• , - ' ~q: 


The re fore , 


+00 

== / l(t) dt 

v/sso -'do 


Applying* condition (AI.S) 
+a; 

/ l(t) dt =1 


(AI.6) 


i.e., the iCiltcr response should be nornolized so that 
condition (..11.6) is satisfied. 

In tliG abBeiice of informtian eibout "lows" . . .• 

H(jw) 

Pron ( Al ,2 ) it follows th.it . under such cond itions , ary 
snail error in gradient reg.oneration' or .noise from any , 
other source causes enornous dis turbunoe , particularly, in 
the low fre quency regions, . . ^e sene e of IPP in the, trans- 
mitter 'helps to reduos such noise generation.. 



Appendix II 

S01.IE QUADIDiViriS E3SUITS 
OOilKSOIINa 

CONTOUR THHBSHOII) AID LOW PASS CUT-OPP 


The follovi/ing analysis although not rigorous and 
r.iay, in fact, contain quite a fev/ loopholes', tries to 
bring out tho close ooiii'.oction botwcon the threshold set 
up for contour extraction and cut-off frequency of the 
filter transmitting tho "Lows". The analysis is done on 

I’igitisod signal of single diraonsion and this can be 

1 ■ ■ 

CiV.i.'i.ly extended' to cases of higher dimensions. 

Lot x(t) be tho time limited sigruil sampled at 
instants t-j # t 2 > • • • j t^, • Then 

x(t) = £ (AII.1) 

i ■ 

where u(t--t^) is the unit step function >/ith discontinuity 
at tj^ and 

« x(t^) - xCt^g^) . 

for all i = 1 , n and xCt^) 0 

Lot x*(t) represent the signal S a| u(t-tj:); 

, . i- 

\7hcrQ .. es'/a^ ,. whenever- 

O' ;'.--other77ise 

i.o,, x*(t) represents the inf prnatipn left out after , 

■tr one nit ting all '-the:, oohtour pp-ints.; ; ;; ; ' 



Tallying the lapliace transforra 


r^(s) 





exp(-stj^) 


(All. 2) 


Substituting s = jw in (AII.2) ws see that 


X»(3w)j 


< 2: “i 


1 exp(-;iwtj^) 


4 f 


(All. 3) 


Ihus the spectrum of the left out portion of the sighaQ. 
has an upper bound given by (All. 3), 


Since the purpose of the "Icy/s” is to supplement 
the information providocl by the contours,: "Lows" should 
contain, as much as pesoiblc information relating to 
x*(t)'. . . 

If w- is tho cut-uff frequency of IPF and e( 3 w) 
represents that portion of X*( 3 w) not contained in "Lows 
cither, 

o(3w) = I*(3w) ; IvYthi^WQ : 

: S3 0 ^ other'-Yise 

lotal power Eq contained in otjw) .is given by ■ 

■ f 0(1!V) iW : 

:| 5 ; f i:X aJ)/V aw;, : ; (AII>4)^ 



IiG't! represent the naxlnuia frequency content of. the 

x(t ) , then 


S 

■^0 


( (% al)/2n ) r 

i 


nax 


dw 

w 




( {X 
i 


c.*)/2^)ln 

i 




pax 

"v; 

0 


(All. 5) 


X a! 

i ^ 


is directly related to the threshold in the sense 


that when tho threshold is increased, a* also increases 
!Proi:i (All. 5) it is evident that an obvious way to counter 


the increase in Bq is to i:icreasc w^. Thereafter, the 
interdependence of these two quantities autonaticiilly 
follov/s. 



Appenclix III 

SYNTHESIS 01 3BC0NSTRUCTI0N FILTER 

Let h(x,y)be tho response of the reconstruction 
filter and b(x,y) the iriage signal. The entire system 
bCiTring tho thresholding operation, is represented by the 
equation 

( T b) * h + b 1 s= b (AIII.I) 

■rfhoi’e * indicates a dot iroduct convolution. The Fourier 
traneforr:! of this equation yields 

Hy) B + BL = B 

where L is tho frequency response of the low pass filtoar. 

For frequencies at which B ^ 0 

^ + jWy = 1-1 

Which When inverse -transformed yields 

7 .h = ujxjy) - l(x,y) (AIII.2) , 

Where :uQ(x,y) is the Dirac delta function in two dimensions. 
If the filter 1 is circularly syEm:etric,CAIII.2) reduces to, 

■the form' ■ ■ ■.■'. 

.. ■V,E(r.) » .u^Xr) - 'iCr)' ,"', ■ (A'lII;# '3) 

The Qorrespohding difforencG equation is given by 



= iij - l(i>J) (AIII.4) 

Where = T i = 3 = 1 , . 

= 0 . otherwise 

If . the LPP is circulfrL’ly symetric , 

Then ocj.n.(AIII.4) bccor'es 

l^(i,j) - h2.(i,3-l) + h^jCj,!) - hjj.(3 ,i--1 ) 

(AIII.5) 

Eq.n. (AIII»4) nust hold goud f or all i anel 3 oven though 
i(i,3) as 0 for 1,3 > n 

But it is desircblfj to hcrwc h^ zero outsiuc the spatial 
extent of liPE so that 

l^(i>3) « 0 i,3 n ' 

It can bo shown that in suoh a situation 
h2,(i,n) ;» 0 ' 

h3j.Cl,T) « (1 (AIII,6^^^^ 

Setting" 

\(i,j ) = 0 ; i ?3+l: ; (A111.7) ; 





A c'otidlocl description of this re construe ti on ^ 

6 

by on Gxanplo^ is 4;; iron in Grahaij's paper . 




Ap..)cnclix rv 

GALCUIu-IIOlT OP B.:.IIDWIDTH ESDUCTIOET EilTIO 


Tho strategy adopted fer conputing the bandwidth 
rc quire me nts of the proposed schece is briefly described 
bolows 

The picture sli-wii in ligure 4.10 is processed 


uoing the pruposod sohoiic 


well as Or chan’ s schene * The 


P.,'.i;4vl('’.th requironontrt ..f both the schonds arc conputed 
u.jing the results ubteihed. Assuidng that 'che rrtio of the 
bandvi dtli roquirononts will not ch^aigc appreciably ■. 
v/liilo processing, nulti-flcvol pictures, the . bandwidt h, re quire - 


nents of the proposed schenc are deduced fron the results 

6 

obtained by Grahnn « The conput ‘t ions leading to the final 
result are indicated below. 


G-rahnn’s Schenej 


ITo, of bits required for - 

(a) Starting point of the contour 
(for 64x64 picture) 

(b) Position inforDatioii. of each 
additional .contour, point; I 

( 0 ) ^G-radiont directioh at each point 

(d) Gradient nsgnitude, for ^each;Cohto^ 




Description of 
fer-tures 


ITuniber Bit 

pre sent re quirenent 


18 18 X 12 = 216 

552 352 X 3 =1056 

6 X 18 = 108 

352 X 3 =1056 
2436 

0}hc Brop.'sod Schenej 

Hoforrins to Table 4, all the approximated contours 
put jo,';othor need only 548 bits. let us see how this 
rit',;uro has boon arrived at . 

OonsiCioring the 5th cntr;>?' of Table 4, 

Uo, of knots =? 4. 

Bits/additional knot = 4 

Ho. of bits for starting point = 12 

end dcrivntive.s= 12 

3 addl. knots = 24 ^ 

, Total Bits - =. ^ ^ ' \ r :■ 

This does not take into consideration the no* of 
bits rcqulrod for code word specifying the bit configuration 
of each cent our. Assuming that each such code .word requires 
6 bits, nn additional 72 bits are required for 12 contours . 
In i'dditioi.'. 36 bits are required for transaitting the^^ ; ■ 

' gj'ad isnt '■ magnit ude , ; 


lie. 

1 Starting point 

2 Additional contour 
point 

3 Gradient magnitude 

4 Gradient direction 

Tot -xl Bits « 



Thtro are 202 ncrc contour points in 25 contours 
ti.'' be transnitted by direct contour coding using the 
proceduro described in Oh..?.pter III, which requires only 
2 bits of each additional node. Ihese points need about 
729 bits for transnission. Another 50 bits are allotted 
for deliciiting code words separating various sections of 
the coder output. Altogether the proposed scherie requires 
548 + 72S + 36 + 72 + 50 = 1435 bits. 


V 


Ihe improve nent 
Gr ..".hari ' a scL’-ene = 


factor of the proposed scheme 
2435/1435 = 1 .7. 


Since Grahan’s sche;x shov/s an inprevenent over 
6 tit POM by a factor varj^ing; from 4 to 23? this scheme 
is likely to give ■; reduction ratio ranging from 6 to 40 
\7lien compared to 6 bit PCM, Prom the results obtained by 
Graham, the low pass signal needs only 1/1 0th of the 
entire bandwidth and hence it can be, safely assumed that 
the same figures hold good for tho entire signal. 

As seen from the results, the approximated contours 
representing well over 400 points need only 656 bits as 
compared to the 729 bits for. 202 remaining points, Ihe 
imp'ortancG of approximation techniques in this application 
is glaringly obvious from the above figures. p , 



o n n nn 


APp"'-'DrX \/ , 
si<":ULat.icij p^:ocr*-Ai';.s 


* IMULAT ION OF LO'-' PA55 FILTER 

.-i.FTc 

■ ' c ■■ . ■ ■ ■ . 

c 3^“ IKTr ar.av-pk 

r Hii , 7 I ‘ Is STCRiiD I|M PIClOW 

!- ..'^AY rT''‘'M'T' AP“R0C5IATi FILTER SAr!PLE 

3 £iy I I' ,.,3v-:3. ■;^'’ of the picture 

C CA>,..;,.:!pi,,.c^. VL, -I, -RPOUEMCY of the filter 

cr-Ex.riop fJluJ” 

I'.TCGFr, PIC.,.piCLn'; ,FIL1 ,FI,lP 

RLAD 3. tPX »■■'¥ 

1 For''AT(?I'») 

READ ? j ( ( P IC( I , j ) , j-3 ,.Mv ) , 1 = 1 ,;:y ) 

2 EcrvAT (A jI?) * ” 

^^?I^'T hi 

111 Fgp;^f!AT(*1 original pICTUPE*/) 

CALL OllTPUT(O) 

OO T 1=1 ,f!x 
DO 3 J=1,MY 

3 PlCLGVd ,J) = o„ 

So, '■•' = ’!„ 

c . 

C FIlTi-R SA^'|-Lt:S— Gtf'EpAT lOi'! THPOUCH THEI;'' '<t-.Ri’jELS 


'll 1 = 1 , in 

R=(I-l)*(I-l) 

'•< = R/( ;>o-!;-SICVA) 

FILL! I ) = 1 »/LXP(R ) 

4 SUf!-.FUM+F ILL( I ) 

SU.V: = 2 --!«-FUf'!-l- 

0051=1,1' 

F I L 1 ( I ) = F 1 1.L ( I ) #1 6 EnO /SUM 
5 FIL?(I )=FIL1(I ) ■ 

12 Cali lp'PF it (nx»ny) 

Do 61 T = 1 »f.'X 
■ 00;61 >1 ,;-:y ' ; ■ ■ 

61 PI CLOl’.' (I , J ) *P I CLOv ( 1 fj ) / ’’ f'.r ,':p n,<-, 
51 FCPMATlX.POlAy 



ouTputt^-g of LC: pass skval 


't' V.'LUTIG i 

nr .r.;. ^ICTjfiu A^:^A,y 


Cali. O! itpuk i ) 

'"\*"siG«r; PICTURE./ 

, I. S', ;= I <;■..-+ 1 

•^UMCH f3 

53 Frr>vAT(2r)^,,4HO'/LrT) 

5 -J STC:-V . 

CfvD 

■’flSFTC LO'APlt 
c Tnl^ 
c 

O'LT IJ' L L,‘)1;.F I |, ( \,y ^|•.;Y ) 

1 :' ; T : ■ G u r - F' ^ ;'■> J C ^ ^ F I [ 1 L ^ ^ ^ ^ '•'’ * *5 0 » 6 5 

i\i=1 - 

■^C ■■' -■ T = ls,V,v 

2 '' j=l ,MY 

■ Nl = I 

1F( j ^.gtuN) ■•:i=fyi 
H2=J 

IF( JuGTcM) .''2 = v 
D o 21 i<r = l ,M 
Il=I-K+l 

addi.=:). 

DO 2.5 ! =1,M? 

Jl=J-L+l 

25 AODl = ADDl+Flui( L )^tPIC( 1 1 , ji , 

■>. vyj. 

I F ( ^'3 0, gF e. F ) N3 = M-i 
IF ( '-'J t.|. FpO ) Go TO 2i 

DO 27 Lrl,N3 ' . 

L ] =I,.H* 1 ■ 

L 2 = J+L 

^^!:!'^I='^DU1+FIH (LI )^fPIC( Ii,L2) 

21 PIClOV ( I ,J)=:plCLOv-'( I,J)+FiL2(i<)*ADDl 

N4»GX-I 

IF { MA oGLoA-' ) .'''AsN-l 

IF(M4«i["</;:) GO TO 2n 
DO 22 K^UWA ■ . . ■ 

. 1C1-='K + 1 ■ 

K?*I+K 
AD'D1 = '-'-' 

;■■ ;Do,V23'; l=i-»^'2;' ■'■ -g g;/ 

■ J l'*-j--L + 1 , 

23 '^DDl = ADfn+Fi:Ll(L)>PIC(K 2 ,Jl ) 



U U U U U U u u 


Ll=L+l 

L2-J+L 

?A ADDIsAl'Dj+F ILI ( Ll ) -Sf-P I C( K2 ? L2 ) 

22 P ICLCW( I ,J) =PICLOW( I ,J)+FIL2 (Kl )*AD[5l 

20 Continue: ' 
returm 

Ei'.iD , . 


■jS' «-****** - sf -K- ^5- jc -X- S -Ji- if-ii- ■«• if- *■){• -if- #•!(• -K ****** •jf 

* 

* Contour extraction ^iooified schEme * 

*' ■ * , ■ 

the Input picture is read into array pic 

THE Contours extracted from the picture arE stopho in x and y 

z Contains the gradient across these edges 

COMMOrVPlCl/PlC( 50»6? ) 

COHMCN/CONTl /MADDOM ( 50»65 ) jX ( 100 ) » Y ( ion) 

Dimension alin(62) 

Integer X9Y»cct,C5T 
CCMMON/PARI/IDIF 

integer pic 

Data BlanksDot/ih ,iHp/. 

READ 1,MX»NY 

1 Format (2 13) 

Do 5 I= 1 »NX 
Do 5 JslsMY 
MADD0M( I+l ,J+1)=0 , 

P I C ( I+l » J+1 ) =0 
5 MADDONd ,J)=0 
PIC( 1,MY+1)=0 
PIC(NX+1»1)=0 
MADD0NM»NY+1)*0 

MADD0N(NX+1»1)=0 4 . 

READ 2, ( (PIC( I ,J) ,J = l,NY)il = J 5NX) 

2 Format (4012) 

I'DIFsl 

call CoNTO(NX»NYd»P) 

32 CONTINUF. 

25 PRINT <9 

9 Format (-sfrl CONTOURS OF . 


nr'i n 


77 


Do 8 I=1»NX 

, DO '6 Jarl jNY , 

ALIN.( J)-E'LANK 

6 IF(PIC( I »J) „£0<.l ) AlIN(J)=DGT 
PRINT in, (alI'N{ J) » J = 1.NY) 

C 10 fORMATdH ,,?0X»55A2) 

10 PORMATdH ,2CX,65A1) 

8 GONtINllE ..... 

. . STOP .' 

End 

SIBFTC COiNTO 
C 

C 1ST. “ COMTOIJ^ STAoT. threshold 
C ICT *- contour CONTI f iUAT I Of-'! THRESHOLD 
C NXjNY -SITE OF THE PICTURE INPUT ARRAY 

SUBROUTINE CONTC{N.X,NY,ICT,IST) 

COM'vON/RIC] /PIC{50,6.5) 

COHfoON /CONTI /HADD.Or-i ( SO ,65 ) sX ( 100 ) ».Y ( 100) 

integer V,Y'.CCT,CsT 

INTLCER ZdC'-) 

IMTEG'EP PIC 
COMHCN /:•■ AP / I DN I » I ONS 
i-X = ] 00 
CCT*ICt 
, CST:=IST 
NXlaNX+1 
NYl-MY+I 

EVALUT 

evaluation of SrADIEi'iT at each POiMT OF THL PICTURE 

Do 1 I=1»NX1 
DO 1 J=1»NY1 
K2 = C ■ 

IF ( J„EQd) GO TO 30 
Kl = ^' 

rF(PlC( I ,J) o£OoPIC( I ,0-1 ) ) C-O TO 31 
K1 = PTC(I ,J)-PICd ,J^1) 

31 IF ( Xl wLT'^CCToAnDoKI cCTo-CCT ) GC TO 32 
lF(Yi..LT-.-.0) K2 = 2 
IF(Kld.T«0) K1*-K1 
KL=MADD0N( I ,J) . 

Call gtgpE(kl,ki ,k2 ) 

MADDON( I »J)''KL 
IFUuEQcNXI) go to 32 
KL*MADDONd-t-l> J) 

K.2sK:2+4^^2'. ■ 

CALL ST0RE(KL»K1»K2) 

MADDf)Nd + l*J)=KL 
' 60 T.0' 32 . ■ . 

30 Kl'spi.Cd,J). ■ ... 

GO'-TO 31 ■ ■ ■ 

32.'K.2«1 ■' 

lF(I,Ecid) GO TO .33, ■ . . 


n r« r; r- ri n 


K 1 = C 

IF(PTC( I » J) »eOoPIC{ I-l ,J) ) GG TO 34 
KX = PIG( 1 ,J)-PIC 

34 IF(Kl9LTuCCToANDoKl„GT6-CCT) GO TO 1 
IF{K1„gEoO) K2=3 
iF.(KloLT.n) Kl = -K 3 
K2S4+K2 
KL=MADD0N( I,J) 

call STGPE(KL»,<1»K2) 

„ !v1ADOON( I ♦ J) =K.L 

IF( JoEOoNYl) GO TO 1' 

KL=mADOON{ I s J+ 1 ) 

K2=8+K2 

CALL store (KL»K1,K.2 ) 

MADDON( I ,J+1 )=KL. 

. Go TO 1 
33 Kl=PIC(r»J) 

GO TO 34 
1 CONTINUE 


EDGE pO*NT * 1 D;'S KEEPS TRAC'< OF THE CONTOUR 

direction at its starti^-’o point 
iL.f.,u Contains the possible dire-.tion of the next 
edge tracing is PcF:FOP’'-^FD By the following section 

Do 20 I = 1 »N.X 1 . 

DO 2 '^ J = 1 

K =0 

INDX =0 

lsw *0 

lDNSs =0 

18 CALL sSTARTd ,J,K,IDnS) 

If ( k e q m 0 w a d a i s w ‘' 4 C:o 1 )G 0 T 021 
IFIKoEq.O) go to 20 

INDX=rNDX+i 
D 1 = I 

L 2 * J 

3 X( INDX)=L 1 
Y( INDx)=L 2 
Z( INOX )=K 
25 K 2 sIDNi *4 
Kl»n 

KL=MADOON(L 1 »L 2 ) 

CALL ST 0 RE(KL,K 1 9 < 2 ) > 

MADDONILl >L 2 )=KL 

Il'Ll 

I 2 *L 2 

IDNT*IDNI +1 

GO TO { 132 » 130 » 133 » 131 ) » IDNl 
134 KL=»MADDON( II ,12 ) 

IDNI = I 0 NT -1 
< 2 »IDNl +2 

lFtK 2 .GE« 4 ) K 2 =K 2 - 4 .: I 



Kl=0 

call ^TGRE(KL*K1,K7) ■ 

MAODO^'^Mo T2)«.<L 
PIC(L3 jL.2 ) = ! . 

PlC{n,I2) = lo 
IF( ISWovIOol) GO TO 6 
.2 11 = 11.." 

, L2=T2 

CALL START(L1,L2»K»IJK) 

26 IF(KocOv.O) GO TO 10 
lNDX=IKiDX+i 

IF( ISW«EQ«1oAMDc InOXoLTvLX) go TG 25 
IFMNDXoLToLX ) GO To 3 
X(INDX)=L1 
Y(INDX)=L2 
Z(INOX)=X 

8 FOR^''AT( 13) 

9 FoRMAT(/6(?.I5,I8»3X)/) 

7 FORMAT (2AI 3) 

PRIMT B,INDX 

PRINT R » ( X ( L) » Y ( L ) 9 ?. ( L ) !. L = 1 » I NDX ) 
PUNCH «i.Ir.:0X 

PUNCH 9 ( X ( L ) ,Y ( L ) 9 Z ( L ) 9L = 1 s IMDX ) 
INDX=1 
ISW = 2 
Go TO 

130 I2=L2+1 
GO TO 13A 

131 I2*L2-1 ■ . 

Go TO 13A 

132 Il=Ll+l 
Go TO 13A 

133 U^Ll-l 

iF(IloFQcC) PRINT 211 9MADDCN(L1»L 
211 format (1 15) 

GO TO 134 

10 IF ( InOx»GT«1 ) GO TO A. 

1 5W=0 

IF{ TDN.‘^«LTc4) GO TO IS 
60 TO 2'' 

4 IF( ISWoEQoO) 60 TO 5 
21 ISWs''' 

PRINT 8,INDX 

PRINT o, (X(L) »Y(L) 9Z(L) *L=1»INDX) 
PUNCH R,INDX 

PUNCH 7» (X(L) vYtL) 9Z{LUL=1 »TNDx) 
INDX=0 

IF{ TDMS‘-LT<.4) go to 18 
GO TO 2C , 

. 5 . rs.iv«i ■ ■ ' 

' Li=i: 

. L2-J , ■ ■■ 

■ Ks-ZM)" 


lO 



GO TO 18' 

6 Do 17 LK=2,INDX 
L<l = T^lr^Y-L<+2 
X(LT-1)=X(LKI-1) 

Y(Li<l)rY(LKl-lJ 
17 Z(Li<l )=Z (Lk'l-l) 

X( 1 ) = Ii ■ 

Y( 1 ) =1.;? , , : , 

2.{1)=-K 
Go TO 7, 

20 CONTIwi!" 

PRINT 11 

11 F0PMAT(^<-0, OVER*) 

■RETURN 

■ ■ '£ND^ 

SIBFTC STORE . : . 

C- ■:■ 

G keeps track, of the nEIGHBoURIN'c edge POINTS 
c THE RElEVENT I.NFORN'AT ION IS CODED By STORE ■ 

subroutine ST0RE(KL,K1,K?) 

K=K2/4 

K3-K2-K^^ 

K = K^-8 

K.3 = K3*?^hS’ 6+K1 
N=(2-«-*K)^(2«-*a) 

1 KL = KL-KL/2-«'*;<*2'’^-«-K+kL/N'«'N-!-K3^!-2^-‘->5K 

return 

End 

SIBFTC START 

C ' ' ' 

C CHECs For an edge point on the contour' 

C IMD - direction WoRct prEVIOU-S EDGE POINT 
c gradient of the neu edge point 
C ■ 

subroutine START( I ,J,K,Ir-'D) 

COMNON /CONTI /MADDON (50 ,6S ) ,X( 100) »Y(100) 
COMHON/PAR/IDNI » IDN.S 

Common /CONT/N x» NY, ccT,csT 
Integer cct,cst 

2 IND1=ImD«-8 
KL=MADD0N( I»J)/2**IND1. 

KLl=KL-KL/256*256 

KL2=KLl/2*^fr6 

KL1 = Kl 1“XL2-?^2**6 
IF(KLUNEuC) GO TO 1 
5 iNDsIND+l 

IF( IND.LT«4) GO TO 2 - 
■ GO', TO 3' ■' ■ 

1 IF{IMDnNE.rL2) KL1=-KL1 
iF(KLH^KuLToO) GO TO 5 . . . 

■ , ,:K=KL1- ', ■■■ 

GO ■ TO .4 ■ ■ ' ' 

3 .■K.*0. . ' " 
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, return 

■ , end'. ■ ■ . 

C ,,•}?■' 

c . , .. CONTOUR. .SMOOTH IMG 

C ¥: 

#IBFTC CONSMO 

C '■■..".■'■'■;■■.■'■■• 

C this ROUTIMt OIVIDES THE CONTOURS INTO SEGMENTS WHEREIN 
C the X*C6CR0ImAtE Is AlVMYS InErEASINGuTHAT is » each 

C : SEGmEwT CAN 6E rEpoESENTEO 3y ASINGlE VALUED FUNCTION 

C X- IMPUT AFRAY CO.NTAINIKG THE GIVEN CONTOUR 
C Y* OUTPUTAFRAY CONTAINING THE SMOOTHED CONTOUR 
C " 

C0MM0M/PARA/L1»L2 vL?. 

COMMON/INO/XOi'lOO) >Y(3!>1C0) 

Integer x,y 

807 READ 800 »LX 

li»L2»l3» Indicators defining tiie starting and Ending 
PORTIONS OF various SEGMENTS. 

READ P.'- ( (X(I»J)s.I = T»3) 9U = 1»LX) 

B . 1 FgpmAT(6 (2I3»I6) ) 

L2 = l 

S-'.O FqpmaT (18) 

LX2=LX-1 
814 LX1=L1+1 

Do 812 I=LX1 »lX 2 
IF (X (1 , 1) -GT«X( 1 , I + l ) ) GO TO S15 
812 IF(X(1,I ) sLToX(l>I + l) ) GO TC 825 
E1=LX 

call SHI ft (1) 

813. Continue 
ilast=i 

ISTARTsI 
L = (^ 

K = 100 

Do B02 1 = 1 , L3 
iFtLoLTol) GO TO 805 
IF(XU ,I )oEQoK) Go to 802 
IF( ILAST-Eq. ISTArT) go To 
IF(L<.Eo.1) go to 805 
IMID=(ILAST+ISTART)/2 
DO 8-^8 J = l,3 

806 Y(J,L)=X{J,IMID) 

805 LsL+1 
Do. 803 J = 1 ,3 

803 YiJ,L).X{J,I )' 

I ST ARTs I 

806 KsX( 1»I) 



u 


802 IlAST=T 


I FM STAR Tc EO«r LAST) GO TO BIT 
IMID= ( ISTAPT+ILA8T)/2 

L “ L + 1 


IF(LoEo„ 1) L=2 
809 J = l,3 

Y(J,l-t)=X(J»IMID) 
809 Y:(J,L)=:X(J,L3) 
all contimue 

PRIMT 800, L 


PRINT 8C4, ( (Y( J,I) ,J=l,3) ,I=l,L) 
804 Format (IH »7(2X,2I4,I8) ) 

PUNCH BOO, L , . 

PUNCHROl , ( f Y{ J, I ) ,J,] ,3) , ,L) 
IF(L1 <.lT«LX- 1 ) GO To 814 
TFtLio-FOaLX-l ) GO To 820 
, Go TO 8''? 

815 NsZ' ' ' 


Go 816 I=LX1,LX2 

816 lF(X(i„noLT;>X(l,I + i)) Go TO 017 
819 L1=LX 

CALL SHIFTCN) 

Go TO 81 ^ 

817 L1=I 


Call shtftcn) 

Go TO 813 


825 N=1 

DO 818 I=LX1,LX2 

818 IF{X(1,I)«GT>.-X(1,I + i) ) Go T:> 

GO TO 819 ■ ‘ 

820 L=1 

Do 310 J=l,3 

81Q Y(j,i)-X(J,lX) 

PRINT Sr'^jL 

PRINT af'4 , ( ( Y ( J , I ) , J=i ,3 ) , I =1 „ L ) 
PUNCH B00,L 

PUNCH 8^1,((Y(J,I),J=1,3),i=i,l} 
GO TO 8^7 


END 

SIBFTC shift 

0 QMMOlT/ino/X ( 3 , 1 00 ) , Y ( 3 , 1 00 ) 

C this routine shifts Each segment to PRE a<;sIGnEu> 
COMMON/PARA/L1»L2,I 3 
integer x,y 
IF( L2oFQ..i ) Go TO 2 
Do 1 I=L2,L1 
IlsI'L?+l 
DO 1 J=l,3 . , 

1 X(J,M)=X( J,I ) 

' 2 ■ L3=L1~L2+1 
: ■ ■L-'2=L1+T 


LOCATIONS 
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IF(N«Eo«1) retupm 

L4=L3/.? 

Do 4 J=l,3 
Do 4 1=1, L4 
N = L3-I41 

2:=x(j,i) 

X{j,n=x(j,N) 

4 X( J,|\i)=Z 
REflloN 
End 

* •»{■**•«•* ** -Ji -»# *#****#-)!.# .Jf. ******** ^{. # ^ 

* * 

* spline approximation problem * 

* * 
«••«*■»■«•*******•*********** *********** 

IBFTC APPROx 

This DETER'^iINES THE mInImum NUMBER OF KNOTS 

needed For a given fit after optimising their 

POoITI.oN ONLY. those ROUTINES WHICH HAVE BEEN , 

Modified have deen listed. « for others refer 

To REF,, 3 AMD 4 

XI L “POSITION OF KNOTS 

ERRL2 “ L2“ EpFOR 

ERRLPR “ MI.NMAX .ERROR 

For other details refer to ref.. 3 and 4 

CoMMON/INPUT/LX,XX( 5C) ,U(50) ,jADD,ADDXIi;.20) ,M0DE 

COMMON/OUTPUT/UERROrI ,FCTL( ‘30) ,XIL( 2? 1 ,COcFL(21 *4) , 
*VorDL<? 292) »KN0T,LMAX,INTE9V 
CoM^'nN /OTHER/LXI ,LXli»'-Xl2909CHANGE,ERROP5ACC5XI (28) 
C0MM0M/C0MS/PHr2 
CoM(^ON/OUT2/ERRL1 ,ErRL2»ERRL99 
COMMON/INTERN/YKSO ) ,¥2(50) ,Y3(50) 

COMMON/ADDN/KX 
Data blank/2h / 
integer Y1,Y2,Y3 
PHI2 = 2c'>^ATAN do ) 

GIVER.=t 0.5 .. 

GO TO 1705' 

1 CONTINUE 

1705 ;ACC=,Oo't,' ■ ■ ■ 

1702 F0R'''>AT(T8»A2) 

. . ;REA.D; 17.02, LX, AMARK . 

■■K'X*.0 

.■VTF'.CLXAT'B.TJ ■,x:.NT . 

. lB(BX*GTo7:oAiND0LX»EB»12) KNTM^ 

1701’' FORMyVT.f2'4l.30 

REAt?<^t,aldYl( I;) ,Y2(I ) ,Y3( D ,I=1,LX) 

■ O0'TO'-23 

1706 F0RMAT<;4:/* 'DATA NOT SUITABLE, FOR APPROXIMATION NOo OF POIf' 

T". 1 .*,I50V'/'.'V 



PUNCH 2702, LX 
2702 FcRHAT(I8) 

PUNCH 170i,(Yi{I),Y2(I)»Y8{I)»I=1»LX) 

GO TO 1705 , 

23 DC 1708 1=1, LX 
XX(I)=Y1(I) 

1703 U( I)=Y2(I) 

WRlTE(f, 5611 ) ( I 5XX( I) ,Ua ) ,I = l9LX) 

ADDXI (1 )=XX(1) 

ADDXI (7)=XX(LX) 

XI(l)=xx(l) 

XI (2)=XX(Ly) 

IF {AmApXoEQ,.. blank.) go TO 12 

15 FCRMAT( I8,/12( 13,3X1) 

RUAD(5,15) LXIs (Y 8(I) ,I=l5LXI) 
lKlN=LXl+3 

MODE=n 
JADCsrliclN-l 
00 16 I=3,JADD 
XI(I-l)=Y2(I*-2) 

16 ADDXI ( I-l )=XI { I-l) 

XI ( JADD) =XX(LX) 

. £'PR0R = FXDKNT<2<. ) 

Go TO 18 
12 IKIN=3 . 

M0DE=0 

JADD=2 

ERR0R=FXDknT{2c ) ' : 
l.f3 Do 9 I»IKIN,KNTMAX . ' . 

■lTER-3 ■ ■ 

' Lx,.I2 = I , 
addxh 1 ) =xx(T-MAx) 

Lxil=Lxl2-l 

lxi=lxti-i 

Do 2J = i '.LXTl ' 

2 IF (XX(LMAX) c^LToXI ( J) ) GO TO 3 
' J=Lxn 

IF ( XX (LMAX ) oFQoXl (Lxl 1 ) ) GO TO ^ 

WRITE(6»700) LMAXsLxII 

700 FORHATI* LHAX out of BOUNDS- ITS VALUE = *,T3,*NOo OF MOTS 
IUCED ^^,13) 

51 IF(ERRL99«GT:»(2o5*GIVER)) GO TO 27 
IKAX=FlOAT(KNTMAX)/1o 5 
IFlLXll -^GT oIKAX) GO TO 7 
■ ' . GO TO 27 

C 3 IF( ABS(XI ( JTrXX (UMAX) > oLTdc-O) GO TO 2.0 
C IF( A6S(XX(LMAX>--XIfU-l> loGEdoO) GO TO 5 ; ^ : 

C 20 IFIABSIXI ( Jl-xrrU-lI ) »bTdoO TO 6 

,3 ADDXT(l)=INT((XTiJ)+Xl(J-lh/2oO) ■ 
lF( ABS{XnU)-Xl {J-1 ) ) «LTo2o ) GO TO ;6 ■ 

5 '.DOrA K»0*LXT,l--." I'l'-' , 

' ' ■■ K1'*LX I2-K:+J 
4 xni<i)»xi(ki-i) 

X I { J) »ADDXi n 5 . 
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!<MOT = J , ' 

85‘J '■CR^iAT ( -If- LMAx >XX { LMAX ) 9 ADDxl ( 1 ) ARE PRINTED BEL0VJ*/?2X » 1 6 » 2^12 
Do S L=J9LXI1 
L1 = L-J+1 . 

a addxi (ij )=xi (L ) 

;'‘'ODF = l' 

JADD=Lyl2-J 

11 Er.RORaFXDK.NTCn, ) 

CALL OwEEP(ITER) ' 

,vrDE«l 

JADD=C 

DiiMB = Fx,PKNT( 2 :> ) 

IF (,EpRLR9-..GToGI VEr), GO TO 9 , 

L>Uf'?=:FxDKNT (-i 6 ) 

IF(£RRLn9,Lt„GIVER) GO To 7 
Go TO' -O ' " 

,6 WPItE(^, 6 ;?.C) I j XI ( J ) 9 X 1 ( J-l ) ^ 

A.?:.' FCRAATf* OPTl''U2ATION DISCONTINUED DUE TO CROWDING OF K'lOTS*// 
iIlE introducing ■k-91^9# -TH knot*//# values of ADJACENT 
2 kNOTS = •':-9.?F12 o ? ) 

GO TO ‘'I 

27 XI (LXl7)=XX(LX) 

LXn=LXl 2 -l 

LXl=LXn-l 

DEL*(XX(L:X)-XX{1 ) )/FL0AT(LXI1) 

Do 26 J=?»LXr2, 

XI ( J-1) = INT.(XI ( J- 2 )+DEL) 

26 ADDXI ( J-2)=XI ( J-1 ) 

MO0E=l . 

■ KNOT ='2. ■ 

jadd=lxi 

DO 28 J=1 jLXI 

28 ADDXIf J)=XI (J+1) 
dumb=FxDk.nt(Oo.) , : ^ 
call SwEcPdTER); . 

MODE-1 

. JADD=0 :. 

. ,DUMB-FxDKNT(2o )d 
9^' CONTINUE' , , 

WRITE (6 975'1) KNTMAX 
7 MODE = l.'".'- 
10 WRITFlEje^'-") 

' jadd=o 

; DUMB = FxDKNT(1o,) , , ■ , . 

LL;',, GO-. TO' 1 '•-AL' .. 

■' A .- T f 0 R'MA T r I 8/1 2. 

, : Sil FpR:MATL^ :GIVEN DATA#//5(IA92F10»5>) 

7,50 Forma T fjfl given accuracy can not be obtained even with ^9:I3 9« 

A^'feL■v■ 'A'''. - 

THIS^ ^ SOLVES the FIXED KNOT . problem : , 

YI 9 Y 2 CONTAIN THE APPROXIMATED CONTOURS 
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For PEtAIlS OW other variable see REF) 
THt OPTION (-1) For arg quantizes the 
UERIVAtES . VORDL{2,2) ^ VORDL ( KNOT , 2 ) 
Cqntaims the End derivatives 


END 


$IBFTC FXDKNt 

Function' fxdknt(arc) 

Double precision trpzv/TjSum 
Logical node? 

Dimension- weight(So) ,cuaerr(50) . • 

COi'vVON />-'Ai:DT/TREND( fO > jTRPZWTCsO) »PRlNT(lOO) 

COMNON / INPUT/LX, XX ('=:0) ?U( 50 ) ■ , JADD,ADDXI ( 20) »MODE 
CCiMHON/0UTPUT/UERROR(5O) ,FCTL(50) 9XIL(22) » COEFL ( 21 » 4 ) , 
■»V0RDL(?2?2) l>KNCT»L^'!AX,INTERV 
COM^jON /BAsIS/FCT(50925) » VORD ( 30 , 23 »2 ) »SC c 25 > , ILAST 
CoMMON/LA3TB/IORDER(28) 9INSIRT(30) ,XKNOT 
COMMOf'70UT?/EPRLl 9EPRL2 9ERRL99 
CommoN/INTEpn/Y1(50) 9Y2{50) 9Y3(50) 

Integer yi!Y29Y? 
coUTVAlFNCE (IpniMT»CHANGE) 

CHANGEsARG 

IFIMoDf .GToO) Go To 29 

xscale=xx(lx)-x>;{i) ■ ■ 

Do JO T = 5 ,30 
1 J I N -S I R T ( I ) = 0 , 

Do 11 l.=l»LX , 

UEPPGRd. )=U (L) 

TRENDf I )-T(XXCL) ) 

11 WEIGHT (L)=W(XX(L) ) 

.. DO 12 L = 2,LX 

12 TRP2WT{L)=(XX{L)-XX{L-1) ) / A ( V/E IGHT { L-J )+WEIGHT(L) ) 

XIl. ( 1 ) =ADDXI { 1 ) ■ 

XIL(2 )=ADDXK2) 

IORD;FR(l)=l 
I0RDER(2)=2 
KNOT=2 
INTERV=1 
Do 19 I-l »4 
■ TLAST.= r 
, CALL NUPAS 
Do 19 L=1»LX 

19 UEPPOP(L) =UERR0R(L)“8CC1 )j^FCT(L»I ) 

' m:odf=i ^ " 

■■ ' ■DO''2P L-i»'LX A : ■ 'D 

20 CUBERR(L)=UERP0R(L) 

7-.F'(.JAP:D = LE-7r)' GO. TO 6;0 
■ ■ DO'-' 2'1-.7=1»7ADD' - - ' 

.21 .ADDXI (. I ).=*ADD.X I ( I,+ 2 5 

-.GO-^^O 51' ■' 

■ 29 GO TO '(4C»40»30)»M0DE 

'".SO.'.'XKNDT'rHANGE' , 

' 'IF'CMODFSO- .go to -' 35 
■ ':-0; MODE-3'=c,TRUE<-.' 



ErSUT1=FXDK'4T 
M0DE=2 
CALL N'UBAS 
KK!0TS\/ = KN0T 
M0DF=3 

GC TO 36. 

35 CALL NUPAS 

36 i-XD'OMT = LR3'jT1-dC( iLAST ) /XSCAlE#BC ( ILAST) 
RETUPW 

AO IF(!<ACt..LT..''-POT5V) GO TO 42 

!<M0T=.i<.M0T3V 

IF ( cNiCTcN:CDE3 ) GO Tq 50 
■ DO 41 L=1 >LX 

41 UEpPOR ( L ) =UERP.OF ( L )-BC ( I LAST ) *FCT ( L s ILAST ) 
Go TO 49 

42 DO 43 L=l9LX 

43 UEPPOP (L ) =CUBERR ( L ) 
lF(KM0TcLEa2) GO TO 48 
iOUM=KNOT+l 

Do 45 I0 = IDUM,KN0TSV/ 

INSERT=INSIRT{ ILAST) 

ILM3 = Ii.AST-3 

Do 44 .K= INSERT » I LM3 

I0RDER(K)=I0RDER(K+1 ) 

44 XlL(K)rXlL(K+l) 

45 IlAST=ILAST>-1. 

Do 47I=5»ILAST , . 

DO 47 L=1»LX 

47 UERROP(L)=UERROR{L)-BC( I )*fCT(L, I) 

GO TO 49 

AS XIL(2)=XIL(ILAST-2) 

IorDER(2)=2 

KN0T=2 

A9 IF(JADDoGT.. >;} GO TO 51 
ILAST=KN0T+2 
INTERVsK.NOT-1 
60 TO 6*^ 

50 IF( JACD.LEoO) GO TO 61 

51 DO 52 10=1, JADD 
XKM0T=aDDXT ( 10) 

CALL NUBAS ' ^ 

D052L“19^X 

■52 UERROR(L) = UE R RO R ( L ) ^BC ( I L AST ) «F C T ( L , I L A ST ) 

60 FxDKNT=DOT(31,2)/XSCALE. 

KMOTSV=KNOT 

61 MODE3=oFALSEo 
lF(IPRlNToLT»0> IPRTNT'-TPRINT 
IF( IPRiNToEQoO) return 
‘GO'tO (70,8n,90),IPRlNT 

70 'WRITE(6,610) 

DO,. 72 1=1, KNOT , 

ILOC=IoPDER( I ) 

'■ D0''.72„ L='1,2'; .■ 

.■ , .'SUMsOaD'^' .. 
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^0 71 J=1,ILAST . 

71 SUM=SUM+BC(J)*V0RD(J,ILGC,L) 

72 V0PDL{ I ,L )=5UM 
Do 75 n'' = l»<NnT 

75 V0RDL(IN>l) = ir.>T(VORDL(IN»l)) 

Al=V 0 RnL(l, 2 ) 

A2 = VOPriL(K.NOT,2) 

LL1 = 2 

NQU=128/LL1-1 

Call QUAiNT(V0RDL(l»?) ,V0RDL(LX,2) ) 

CALL IwTEpp 
call EVAL 
Do 76, IN = 1»LX 

76 UERF;GR(IP)=FL0AT(YL(IN))-U(IN) 

IF(ARGoLToC) go T 0.P,5 

. FxD'<NT = DOT{.31.»2)/X3CA.lE 

fcrmat{ // iH of quahtisation levels 

1 . ,= •?!-jP 12.69 15XFl2o6 ) 

DO 73 t = 1<,!<NGT ■ - 

73 PPR''! 620s.I ,XIL ( I ) sVORDLf I ) 

,610 FoPAAT(42X:5H!<NCTS927X.* VALLES '*//) 

620 FQpf.-AT{^5X..'^Hxl ( , I293 H)= , Fi 2 A, 22 X j F] 2 <>5 ) 

PRTfT 63 0 »vorDL(1,2),VGPOL(K.MOT9 2 ) 

630 FoPMATc^^O*, 35X»* END DERIVATIVES*,* = 2 « 6 , 15X »F12 of ) 

■ 80 Erol 2 = SOPKFXDKMT) 

ERRLI = Co ■ 

ERRL99=n« 

Do 82 L=l9LX 

OlF =ABS(UERROR(L)*k'ETGHT(L) ) • 

IF{ERRL99 oGToDIF) Go to 81 
LMAX=L 
ERRL9q=DIF 
SI Erpli=FPRL]+DTF 

if( arGoLToO) return 
82 Continue 

ErrlI=EFRL] /float (LX) 

Go TO ( 909 <^ 6 s.Q 6 ) 9 IPPINT 
90 Continue 

PRINT i: 09 LNAX 9 ERPL- 9 . 9 EPRLI 9 ERRL 2 
100 FoRMATdH 9*LMAX»ERRL99 = *953X,I59F12:::6/ 

1*0ERRL1<,ERPL2 =. *.d'^..K.?(Fl 2o6dOX)) 

PRINT 8''C'9 {XX(L ) d'CTLCL) 9 L-I 9 LX) 

p.bO Format (*ofinal values of xx^fctl */( 2X9.9 Fi 6 o 8) ) 

V0PDL(KN0T(1»2)=A1 
V0RDLtKN0T,2 )=A2 : 

96 RETURN 

E'ND ■ . ■ 

THIS ROUTINE QUAMTI 7 .ES THE AniGlE O.F . 

Inclination to the required fineness e 
N-.NUmBEROF,.LEvELS '.' . 

A,B r. quantities TO BE GUANTISED 


SIBFTC QUANT 



5uBpoUtInE OUAkT(A,B,n) 

C0M^'0^:/C0^!S/PHI2. 

AN=N , 

, A,\!=:AM/pHI2 
AN=AM/9, 

Al5:ATAW(A) 

I F { ABS { A 1 ) oGTo ( PHI2*2 o/3 o ) ) AN=4o*AN 

Ki=AES(A1*AM)+0,F 

A2=FLaAT (K1 ) /AN 

IF{AB.S{ A2-PHI2ULTcCoOA2) A2 = PHI2-0c-002 
AI-LToO. ) A2=-A2 
A=TAN(A2) 

AN = K' 

Am=AN/pHI2 
AN=AN/2o 
BisATAn (B ) 

IF( ABS(5l ) cGTo (PHI2’«-2v./2c ) ) AM=4 -,-s«-AN 

K2=ABS{B1*AN)+0c5 

B2=FL0AT (K2 ) /AN 

IF( ABs(B 2-PKI2) -LToOoOC?) B2=PH 1 2-0 ^ 002 
' I,F(B loLToO. ) B2 = -B2 
B=TAN(B2 ) 

C * AA -x- -if -X- > V -5^ ^ ^ -:v *?f ^ if -if if if if if if 'X- if if if X- if if if if 

C -X- if 

C if PICTUPE pFGONSTPUCTICM . if 

C if if • 

(” if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if i^ if i** if ^ if if if if if . 

SIBFTC 

COiyMOM/MAlN/PlC( 66 »65 ) 

CoK'-^ON/F lLT/PlCh(65r653»Fl3i <23 ) »F IL2 ( 26 ) ^HX ( 12 » 1 2 ) 
Co^<'^’ON/PAR/NX ».NY 
CoMf''ON/OUTsF/ALlN{65) ?!GR<2) 

Integer pic9pich,fili»fil2 9 1"^x 
Dimension fill(25) 

Integer x(ioo)»Y(iO'^),z(ioo) 

Integer pictn.t(66»63) 

NX =41 
NY =46 
NslO 
FELE=0 
SUM=0, 

SIGMA=4oO 
NM=NX 

IF (NMolToNY) NM=NY 
Do 152 I=1»MM 
DO 152 J=1»NY 
PICH( I ,J )=0 
PICINTT i »J)=n 
152 PiCd ,J)=:' 

;C , ^ 

C H p. F . response. GENERATIOfi ., 

.C ' ■ 

Do 4 1 = 1 »,N ., . , 

R=( i-i):^( i-i ) 

R*R/( Pa-^-SIGMA) 



FILL( I )=lo/EXD{R) 

4 S!JM = SUm+FIlL Cl ) 

SUM = 2o-SK?UM-i, 

00 5 I=1,N 

Fili ( r) = FlLL{ n-:5-l ^tr.3/s‘J''i 

5 FIL2 ( I )=FILi ( I ) 

Do 505 

505 HX(N,I)=0 

DO 501 0=2,^ 

.i = N-J+2 

IF( I GO TO 504 

Il=I+l 

Do 502 Ll-JlsM 
502 HX(I-I,L1)=' 

5''4 Hx ( I-l , r ) = (FIlI ( I )*FrL 2 ( n-i- 2 *HX( < 59 1 ) ) /.2 

>-iX ( I-l ,1-1 ) = H ILI ( I )^5FIL2( I-D+Hx ( I ? I-l )+HX{ I-l »!) 

if.cIoEqo?) go tc '.C?, 

■ 1 2 = 1-2 

-'0 5vl L = 1‘.. 12 
' ^= 1 - 1 - 1 ■ 

5':i Hx ( I-l .K ) =FILI( I I'K-r'IL?. (lO+HX ( I ,Ki: 

5;‘3 HX{1,I) = (1'^ ’■'f'OO-F ILI { D^FILZC 1) 8/4 

75 FORMAT f T 8 ) . ' 

80 READ 75<.LX 

IF( LXr-EOo' ) GO TO -^8. 

READ 76 9 ( X( I ) ,Y{ I ) sT; I ) »I = 1‘. LX) , 

76 FORMAT(6(2T39i6) ) 

, Do 77 T = 2 ,l.X 

Kl=X( I } 

K?. = Y( I) 

K=X(. 1-7 ) 

L=Y ( I-t ) 

Kl = X( I-l )-X ( I ) 

K2=Y( I-1)-Y( I ) 

KL=7( I-l) 

IF { I AEs( K1 ) o&To 1 oCRo TABSC l<2 ) ^GToI ) GO TO 2^ 

IF ( IABs(Kl^K2 ) oEGvl ) GO TO 23 
IF(KloEQ.O). GO TO- 21. 

IF (51° EQ 1 “ /'ND n K 2 = o 0 ) L= L— 1 
IFCKI^EO^ (-1 ) pAN'DoKE^EOn) KL=-KL 
. PIC{K,L)=KL' 

GO TO 77 ■ , 

21 IF(K 2 oEQ.- (- 1) ) GO To 22 

K.= K-1 

KL=-KL 

22 PICINT(K,L)=KL 

77 CONTINUE 
GO TO 8^ 

78 CONTINUE 
: InD1=1 - 

85 CONTINUE 
. CALL HIFI-L 
DO 92 . 1 = 1, NM... ■ . 



n o n o 


00 92 

^>2 PIC{ I ,J)=:-^ 

IF( InDi..£Qc 2 ) GO TO 90 
Do 83 ’l=i s.N'P-1 

1 F ( I a Eq..- NM ) GOTO 8'^i- 
Jl=I+l 

Do 83 J=Jl9:NM 

K = PICH(I sJ) ■ • . 

PICH( I , J.)=PICH( J, I y 
PICH{J,I)=K 
K=P,ICImT(T ,j:) 

?ICINT(IfJ)=PICLMT{J,I) 

83 PICINT CJ , I )=K 

84 lMO] = rN!Dl + l , 

Do. ?5 1 = 1 

Do ?.p J = 1 , , 

PIC{ I,j)=C 

25 PIC(r,j,).=PIClMT('j5l.} 

MY = Mx 

GO TO 85 
90 NX=MY 
NY = FM 

READ 7o s ( ( PI C ( I , J ) f J = i ,NY ) » 1 = 1 sMX ) 
c 79 Format (401 2 ) tagore , 

79 Format (401 2 ) 

Do 86 1=1, NM 
Do S6 J=1,NY 

PICH( J ,I )=:PICH( J , I ) /locoroc 
86 PIC(I,J)=PIC(I,J)+PICH(J,I) 

Do 251 1=1 ,6 

• PRINT ,161 

161 Format ( iHi////// ) 

CALL OUTPUT! 3, ) 

PRINT 3.62 _ . 

162 FQp,r4AT(////30X»-K- fig 4011 ptCRiVED pICTURp - GRAHAM 
251 Continue 

PUNCH 79» ( (PIC( I »J) ,J=:1,NY) ,I = 1,NX) 

STOP 

End 

vHs -iv-s ■?<■-)(■ I 

^ HIGH PASS FILTER ' V 

',.' 0 '' .. , 

■! 1 - 4 4Hi- -S- « r- ^ -if ^ ^ ^ 

SIBFTC HIFTL ' ■ - ^ 

subpoutine hifil 

CoMMON/MAIN/PIC (65,65) 

CoMMON/FiLT/PiCHi (65,65 ) ,F IL 1 ( 25 ) ,F IL2 { 23 ) ,HX§[;i|l2.) i 
G0MM0N/PAR/NX,NY,N 

C0MM0N/FILT2/FElE,FC0L (65 ) ,FrcW(65 ) 

Integer gradx, PI CHI, FI LI, fi,l2,hx 
. Integer fele,,.fcol,frow 



Do 3'> I = X,NX ^ 

Do 30 J=1 ,NY P 

N3, = I . ■/ 

IF( I0GT3N) N1=N 
N2=J ■ , 

IF ( j„GTcN) N2=N 

Do ,43 r.=isNi 
Do 35 . L = 1,N2 

D 1 = J-~L 

P I CHl ( I , J ) =:p I Chi ( I s j ) +HX ( K: ? L ) I C { 1 1+T , Jl+i ) 

34 ;\13=:>.!Y-J 

IF(N!3oGEoN) M3=;W_T[ 

IF ( N3 oI_E «0 ; Go To 43 
Do 37 |.,. = 1,M3 
■ L1=L+1 
L2=J+L-1 

37 PIChi ( I ,j)=:pICHi ( I sj j+HX{ KsLl }^PlCMi + i,L2 + l3 

43 ComtINUE 
N4=NX-I+1 
IF{M4<-GT«N) N4=M 
IF< Ki/ 4 o|_£oC, ) Go To 3i^ 

Do 2 8 k: = l »M4 
Kl = K 

K2 = I+K,^1 
Do 39 L=l!'W2 
J1=J-L 

39 P I CHl ( I 9 J ) = p I CHi ( I , j ) -HX ( K1 3 L ) *P IG ( K2 + '! s JT -fT ) 

42 IF(M3oLEeO) GO TO 38 ' ' * 

Do 40 L=1»N3 
L1=L + 1 , 

L2 = L+J-1 

40 PICHl { I , J)=PICH1 ( I,J)-HX(!<l9Ll)^PIC(K2+l,L? + -| ) 

38 Continue . , " 

30 continue 

RETURN 

■■ end ££■':' T 

23 PRINT .4»K1 vK2»K 

4 Format {•«• points not co.nsecutive*/* ki 9^'2 »i< = 

STOP 


END. -OF .SIMULATION: ^ .01 Bjw tWAL. ':y»ig 4yiir^ 


E.£ ^ -S uTi 




